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ABSTRACT 

We calculate the pressure in QCD with two and three light quarks on a lattice of 
size 16^ X 4 using tree level improved gauge and fermion actions. We argue that 
for temperatures T > 2Tc systematic effects due to the finite lattice cut-off and 
non-vanishing quark masses are below 15% in this calculation and give an estimate 
for the continuum extrapolated pressure in QCD with massless quarks. We find 
that the flavour dependence of the pressure is dominated by that of the Stefan- 
Boltzmann constant. Furthermore we perform a calculation of the pressure using 2 
light {mu,d/T = 0.4) and one heavier quark {rris/T = 1). In this case the pressure is 
reduced relative to that of three flavour QCD. This effect is stronger than expected 
from the mass dependence of an ideal Fermi gas. 



1 Introduction 



One of the central goals in studies of the QCD thermodynamics on the lattice is 
the calculation of the equation of state for QCD with a realistic mass spectrum. 
Controlling the influence of a heavier strange quark on the QCD phase transition 
as well as understanding its contribution to bulk thermodynamic observables, e.g. 
the pressure or energy density, is of fundamental importance for the analysis and 
interpretation of heavy ion experiments which look for signals from the quark-gluon 
plasma phase of QCD. While the former problem is closely related to the chiral 
structure of QCD and requires numerical computations with light quarks close to the 
chiral limit |T| the latter question can be addressed in computations with moderate 
quark masses already. 

An increase in the relative abundance of strange particles at high temperature 
and density is being discussed gnature for the formation of a quark-gluon 

plasma [Q. A first analysis of the contribution of the strange quark sector to the 
energy density has been performed already some time ago in a lattice calculation 
1^. This calculation suggested that even at temperatures a few times the critical 
temperature the strange quark contribution to the overall energy density is strongly 
suppressed. The energy density in the strange quark sector was found to reach only 
half the value of a non-interacting Fermi gas. This conclusion, however, has been 
drawn by separately analyzing operators which in a non-interacting gas would de- 
scribe the contribution of light and heavy quarks, respectively, to the overall energy 
density. Such a separation of different contributions is questionable at temperatures 
a few times T^. In fact, the experience gained from lattice calculations of the equa- 
tion of state in the pure gauge sector Q], the failure of a perturbative description 
of its high temperature behaviour as well as the success of hard thermal loop 
resummed perturbative calculations ^ suggest that the high temperature phase 
of QCD remains non-perturbative even at temperatures several times T^. This pro- 
hibits an isolated analysis of the contribution of different parton sectors to bulk 
thermodynamic observables, e.g. the free energy density /(T) which in the thermo- 
dynamic limit yields the pressure, p{T) = —f{T). One rather has to analyze the 
variation of /(T) with the number of flavours and its quark mass dependence to 
deduce the effect of a non-vanishing strange quark mass on the thermodynamics. 

The pioneering calculation of the strange quark contribution to the QCD equa- 
tion of state uses the standard Wilson gauge and staggered fermion actions. It 
still is strongly influenced by cut-off effects and moreover, makes use of partly per- 
turbative relations in the calculation of thermodynamic quantities. Both problems 
can be handled now much better. The use of the integral method allows an en- 
tirely non-perturbative calculation of thermodynamic observables. Cut-off effects 
can be strongly reduced in calculations with improved gauge |P, 0, |lT], |12| and 



fermion ||13| , |14| actions. 
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We will present here results from a calculation with improved gauge and improved 
staggered fermion actions for two and three quarks of equal mass as well as two light 
and a heavier strange quark. Based on an analysis of the remaining cut-off effects 
and the quark mass dependence we will give an estimate for the free energy density 
(~ pressure) of QCD with massless quarks at temperatures T > 2Tc. In addition 
we discuss the contribution of the heavier strange quark to the free energy density 
in this temperature interval. 

This paper is organized as follows. In the next Section we will describe the lattice 
action used for our calculations and discuss the cut-off dependence of thermodynamic 
observables in the infinite temperature, ideal gas limit. In Section 3 we present 
details of our numerical calculation and give the basic numerical data entering our 
analysis of the pressure which is performed in Section 4. In Section 5 we discuss the 
extrapolation of our results to the continuum limit for which we give a first estimate. 
Our conclusions are given in Section 6. 



2 Improved action, Cut-ofF dependence 

The equation of state for QCD with two light quarks has been analyzed recently on 



lattices with temporal extent A^,- = 4 [|T5l and 6 ||T6|. In these calculations, which 



have been performed with the standard Wilson gauge and staggered fermion actions, 
a sizeable cut-off dependence has been observed. The general pattern, a strong over- 
shooting of the continuum ideal gas limit, is in accordance with the known cut-off 
effects for an ideal quark-gluon gas calculated with these actions on lattices with 
finite temporal extent. Also the first thermodynamic studies performed for four 



flavour QCD with an improved fermion action [|T3| indicate that the qualitative 
features of the cut-off dependence closely follow the pattern seen for an ideal gas. 
Analyzing cut-off effects in the ideal gas limit thus provides useful guidance for 
selecting an improved action for thermodynamic calculations. 

In our calculations we use a tree level, Symanzik improved gauge action, which 
in addition to the standard Wilson plaquette term also includes the planar 6-link 
Wilson loop, and a staggered fermion action with 1-link and bended 3-link terms. 



Z{T,V) 

Sg 



I n df/.,/.e-^^« n ( / n e- 



1/4 



(2.1) 



C4 Spidquette ~l~ Spi^nar 

- E \ fl-^ReTr 



3 





(2.2) 



^ J/l Ixl 11 




y X y 





(2.3) 



Here we have made explicit the dependence of the action on different quark flavours, 
/, and the corresponding bare quark masses and give an intuitive graphical 
representation of the action. With ri^{x) the staggered fermion phase factors are 
denoted. Further details on the definition of the action are given in |]I4|. The tree 
level coefficients cf and appearing in Sp have been fixed by demanding rotational 
invariance of the free quark propagator at 0{p'^) ("p4-action"). Moreover, the 1-link 



term of this action has been modified by introducing "fat" links [17] with a weight 
uj = 0.2. In the infinite temperature limit the fat link term does not contribute 



to thermodynamic observables nor to their cut-off dependence. Also at 0{g^ 
effect has been found to be small 



its 
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It thus is expected to be of httle importance 
for our current analysis, which is focused on the high temperature behaviour of the 
pressure. The fat links are, however, known to improve the flavour symmetry of the 
staggered fermion action [T^ . Their contribution thus should become of significance 
in thermodynamic calculations with light quarks close to Tc. 

Even with our improved action cut-off effects are still quadratic in the lattice 
spacing a. In thermodynamic calculations this translates into a quadratic depen- 
dence on the finite temporal extent of the lattice, (aT)^ = Compared to 
the standard gauge and staggered fermion actions these contributions are, however, 
drastically reduced in magnitude. For j3 oo the p4-action yields a much more 
rapid approach to the continuum ideal gas limit and deviates little from it already 
on lattices with rather small temporal extent N^-. The deviations are smaller than 
5% already on lattices with temporal extent Nr = 6. Even with an improved gauge 
sector this accuracy is reached with the standard staggered fermion action only for 
Nr > 16. The cut-off distortion at some small values of N^- is given in Table |l|. 
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PG{Nr)/PSB,G 


PF{Nr)/psB,F 


/ ( AT \ 


^PF[Nr))/T^ 




gauge action 


p4-action 


Naik-action 


p4-action 


Naik-action 


4 


0.9284 


0.5932 


0.8070 


1.3393 


1.5991 


6 


0.9925 


0.9378 


0.7718 


-0.0867 


0.8079 


8 


0.9983 


0.9778 


0.9303 


0.0236 


0.1156 


10 


0.9994 


0.9889 


0.9809 


0.0155 


0.0202 



Table 1: Cut-off dependence of tlie pressure on lattices with temporal extent A^^. We 
separately give the gluonic contribution and that for massless fermions normalized 
to the corresponding continuum ideal gas values. The second column gives the result 
for the tree level improved gauge action defined in Eq. (2.2) and the third column 
is for the p4-action defined in Eq. (2.3). In the third column we give for comparison 
the results for the 0{a?) improved Naik action |TE]. The last two columns show the 



violation of a basic thermodynamic identity for the ideal Fermi gas (one flavour) 
due to finite cut-off effects. 



In the last two columns of Table |1| we also show {tp — 3pf)/T'^ calculated by 
using the lattice versions of standard thermodynamic relations, i.e. p/T = V^^ InZ 
and e = —V~^dlnZ/d{l/T). The resulting integrals which have been evaluated 
numerically are given for a free massless fermion gas, 
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where the zero temperature contributions to Pf/T'^ and ep/T'^ have been subtracted. 
Here the function uj'^{p) = 4X]jl=i PiPp.) is introduced where in the case of the Naik 
and p4 action f{pp) is given by 



f{P,) 



9 1 

— sin(p^) - — sin(3p^) (Naik - action) 
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/(Pm) = g sin(p^) + 2 sin(p^) ^ cos(2p^) (p4 - action). (2.7) 

In the temporal direction only the discrete Matsubara modes po = (2^?.o + l)7r/A^T 
with no = 0, 1, {N^ — 1) contribute. The deviations from zero in {ep — ?>Pf)/T'^ 
are a direct measure of the violation of basic thermodynamic identities, valid for an 
ideal gas, due to finite cut-off effects. This also supports our preference for using 
the p4-action in thermodynamic calculations rather than the Naik-action. 



3 The numerical calculation 



Our calculations for two and three flavour QCD have been performed with quarks 
of mass mu^d;L = 0.1 on lattices of size 16^ x 4, i.e. niu^/T = 0.4. In addition we 
perform calculations with two light quarks of mass mu,d;L = 0.1 and a heavier quark 
of mass rris-L = 0.25, i.e. nis/T = 1.0. To normalize the pressure and also in order to 
extract a physical temperature scale additional zero temperature calculations have 
been performed on symmetric 16'^ lattices. We have used the standard Hybrid R 
algorithm [|I9] with a step size Ar< mu^d^L/'^ and a trajectory length r = 0.8 to 
update the gauge and fermion fields. On the 16'^ x 4 lattice we collected 2000 to 
3000 trajectories for values of the gauge coupling near the critical point and about 
1000 trajectories for values away from it. In the zero temperature simulations up to 
800 trajectories were generated to obtain a statistical error comparable to the finite 
temperature calculations. 

In a first step we determine the transition region for our choice of parameters. A 
pseudo-critical temperature has been determined on the 16^x4 lattice by locating the 
peak position of the susceptibility of the Polyakov-loop and the chiral condensate, 
respectively. The location of both peaks has been found to coincide within errors. 
The signal in the chiral susceptibilities is, however, far more pronounced than in 
the Polyakov loop susceptibility. The resulting pseudo-critical couplings are given 
in Table |. 

In order to determine a physical scale for our finite temperature calculations we 
have calculated the heavy quark potential and also performed spectrum calculations 
at T = 0, i.e. on the 16^ lattice. The heavy quark potential has been determined 



in the usual way from smeared Wilson loops (for details see |rT|). From its long 
distance behaviour we extract the string tension which then yields Td \fo and also 
defines a temperature scale for our calculations of the pressure. The resulting critical 
parameters are also given in Table |^. 

We note that Td \/a shows little dependence on the number of flavours, al- 
though the present analysis certainly is not yet indicative for the behaviour in the 



6 



flavour content 


/3c 






rripsa 




2 


3.646 (4) 


0.271 (10) 


0.480 (10) 


0.958 (2) 


1.377 (25) 


2+1 


3.543 (2) 


0.271 (11) 


0.480 (10) 


0.962 (3) 


1.343 (20) 


3 


3.475 (2) 


0.283 (11) 


0.470 (9) 


0.967 (1) 


1.415 (15) 



Table 2: Pseudo-critical couplings, string tensions calculated at these couplings and 
the resulting pseudo-critical temperatures for = 2 and 3 as well as QCD with two 
light and a heavier strange quark. In addition the masses for the light pseudo-scalar 
and vector mesons are given. 



chiral limit where Tc/ a/o" seems to be about 20% smaller than our current result 
20[| . Furthermore we have performed a spectrum calculation on a 16*^ lattice at the 



pseudo-critical couplings. For the ratio of pseudo-scalar and vector meson masses 
we find in all three cases, mps/m^ ~ 0.7 (see Table This also indicates that the 
quark masses used in our current analysis are certainly too large to investigate in 
more detail the temperature interval close to Tc. In the high temperature phase, 
however, the dependence on the bare quark masses is strongly reduced [|T^ and 
a reliable estimate of the pressure in the chiral limit becomes possible. 



4 The pressure 

The free energy density, / = —TV~^ In Z, which in the thermodynamic limit directly 
yields the pressure, p = — /, can be calculated using the integral method p. As 
the logarithm of the partition function is not directly accessible in a Monte Carlo 
calculation one first differentiates the partition function, Eq. (2.1), with respect to 
the coupling /5 = G/g"^ and subsequently integrates the resulting expectation values 
of the gauge action. 



rp4 



= -(w) / df^'HSGh-iScW) . (4.1) 



Here No- and denote the spatial and temporal extent of the finite temperature lat- 
tice, respectively. In the limit — * oo, Eq. ( |41| ) gives the pressure, p(T) = —f{T). 
Experience from earlier calculations shows that aside from a small temperature in- 
terval around Tc this limit is well approximated by N^j > ANr- Moreover, {Sg)t is 
the expectation value of the gluonic part of the action at finite temperature. The 
zero temperature contribution {Sc)o calculated on the 16"^ lattice is subtracted to 



normalize the pressure to zero at T = 0. Strictly speaking Eq. ([4.1|) gives the differ- 



ence between the ratios f{T)/T^ calculated at two different temperatures, T = T{[3) 
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and To = T{(3q). We choose the latter such that f(TQ)/TQ ^ 0. Within our numer- 
ical accuracy this is the case for Tq ^ 0.6Tc. Calculations of the differences of the 
action expectation values have then been performed at about 15 different values of 
the coupling which cover the temperature range 0.6 Tc < T < 4 T^. 




Figure 1: String tension calculated on 16^ lattices (a) and action differences from 
16^ and 16^ x 4 lattices (b) for nj = 2, 3 and (2+1). 

The basic input for our analysis of the QCD thermodynamics thus consists of a 
zero temperature calculation of the string tension, which is used to fix the tempera- 
ture scale, and a calculation of action differences. The latter can then be integrated 
analytically to give the pressure according to Eq. ( [4.1[ ). These data are shown in 
Figure |1|. 

For the interpolation of the string tension data we use a renormalization group 
inspired ansatz ||21||, 



vW(/3) = R{(3){1 + C2a\(3) + Cia\(3))/co (4.2) 

with d = R{(3)/R{(3). The interpolation parameters are given in Table |^. Although 
this ansatz is, in principle, suitable for an extrapolation to the continuum limit we 
stress that it is used here only as an interpolation for the string tension data and as 
such is valid only in the interval indicated in Table ^. 

The systematic increase in the action differences with increasing number of 
flavours visible in Figure ^ leads to the increase of the pressure with increasing 
number of the degrees of freedom, which is apparent from Figure where we show 
the pressure for Uf = 2 and 3 as well as the (2+l)-fiavour case. In fact, in the 
case of the simulations with two and three light quarks, respectively, we observe 
that this flavour dependence can almost completely be attributed to that of an ideal 
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flavour content 






Co 


C2 


C4 


2 


[3.6,4.4] 


3.70 


0.0570 (35) 


0.669 (208) 


-0.0822 (1088) 


2+1 


[3.5,4.4] 


3.60 


0.0526 (32) 


1.026 (224) 


-0.1964 (1065) 


3 


[3.4,4.2] 


3.50 


0.0448 (15) 


0.507 (115) 


-0.0071 (677) 



Table 3: Fit parameters used for the interpolation of string tension data, 
quark-gluon gas. 
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Here Qf counts the effective number of degrees of freedom of a massive Fermi gas. 
For a massless gas we have, of course, gf = nj. In general we define 



9f= aimf/T) 

f=u,d,.. 



(4.4) 



with 



360 r°° I 

g{m/T) = — / dxxJx^ - {m/Tf ln(l + e"") 



Itx'^ Jm/T 



(4.5) 



For the quark mass values used here one gets (7(0.4) = 0.9672 and g{l) = 0.8275, 
respectively. The correspondingly normalized curves are given in Figure |]b. This 
indicates that in the presence of a heavier quark the deviations of the pressure 
from the ideal gas value is larger than in the massless limit. This is in qualitative 
agreement with the observations made in B. 



5 Continuum Limit: An Estimate 

Our current analysis is restricted to a single temporal lattice size, i.e. = 4. We are 
thus not yet in the position of performing a complete extrapolation to the continuum 
limit. With our improved action finite cut-off effects are, however, strongly reduced 
at high temperature. We thus may attempt to give an estimate for the continuum 
equation of state for massless QCD with two or three quark flavours at temperatures 
not too close to Tc, e.g. T > 2Tc. 
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Figure 2: The pressure for n/ = 2, 2+1 and 3 calculated with the p4-action (a) 
and the normalized values p/psB (b).The arrows indicate the continuum ideal gas 
limits for two and three flavour QCD with quarks of mass m/T = 0.4 as well as the 
case of two flavour QCD with m/T = 0.4 and an additional heavier quark of mass 
mjT = 1. 



The influence of a non-zero, yet small to moderately large quark mass is small 



at high temperatures [T^, |15| and, moreover, seems to be well described by that 
of an ideal Fermi gas. Figure ^d. The ratio of Stefan-Boltzmann factors for QCD 
with two and three light quarks of mass m/T = 0.4 and massless QCD is 0.981 and 
0.978, respectively. Thus, this seems to be a minor source for systematic deviations 
from the massless continuum limit. The main source for systematic errors clearly 
still are finite cut-off effects. 

The analyzes of finite cut-off effects in the pure gauge theory have shown that at 
temperatures T ~ (2 — 4) Tc the ideal gas calculations correctly describe qualitative 
features of the cut-off dependent terms. However, they overestimate their influence 
by roughly a factor 2. If this carries over to calculations with light quarks, which 
similar to thermal gluons also acquire a thermal mass of 0{g{T)T) we may expect 
that the finite cut-off distortion in our numerical calculations is also reduced by a 
similar factor. From Table |l] we find that in the ideal gas limit our improved action 
leads to results which are 26% and 29% below the continuum value for = 2 and 
3, respectively. Combined with the small systematic errors resulting from the use 
of non-zero quark masses we thus expect that the continuum equation of state for 
massless QCD at temperatures T > 2Tc is about 15% above the values currently 
obtained in our analysis. This estimate for the continuum limit is shown for two- 
flavour QCD in Figure |^ where we also show results from a calculation with the 
standard Wilson gauge and staggered fermion action on lattices with temporal extent 



A^T- = 4 and 6 [0. These latter data lie substantially higher which is in accordance 



with the larger cut-off effects for the unimproved actions. 
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Figure 3: The pressure for nf = 2. Shown are results obtained with the p4-action on 
lattices with temporal extent A^^ = 4 (line) as well as with the standard staggered 
fermion action on A^^ = 4 (squares) and 6 (circles, triangles) lattices. Also shown 
is an estimate of the continuum equation of state for massless QCD (dashed band), 
based on the assumption that the systematic error of the current analysis is (15±5)%. 



The lattice studies of the QCD free energy density, or equivalently the pressure, 
indicate that deviations from the infinite temperature, ideal gas limit are about 
(15-20)% in the temperature interval 2Tc < T < 4Tc. This is quite similar 
to what has been found in the pure gauge sector. We note that also the HTL 
resummed perturbative calculations ^ suggest similar deviations from ideal 
gas behaviour. This does form a basis for more phenomenological quasi-particle 
models for the QCD equation of state p3, p5 . 



6 Conclusions 



We have calculated the pressure in the high temperature phase of QCD using im- 
proved gauge and fermion actions. Our analysis focuses on the flavour dependence 
of the pressure and its dependence on a heavier (strange) quark mass in the high 
temperature plasma phase. We find that the quark mass dependence closely follows 
the pattern expected from the analysis of an ideal Fermi gas. We observe, however, 
a significant reduction of the contribution of strange quarks relative to that in an 
ideal gas. Interactions, which at temperatures a few times lead to a strong re- 
duction of the pressure in QCD relative to that of an ideal gas thus also show a 
significant quark mass dependence. We note, however, that the current analysis has 
been performed with a fixed ratio m^/T rather than a fixed heavy quark mass. The 
latter will be needed to arrive at quantitative predictions for the plasma phase of 
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QCD. 



So far we only can give an estimate for the continuum extrapolated pressure in the 
high temperature phase of QCD. To complete this analysis and in particular in order 
to analyze the transition region and the nature of the transition in (2+1) flavour 
QCD further calculations on lattices with larger temporal extent and smaller quark 
masses are needed. As the p4-action strongly reduces the cut-off dependence it pro- 
vides a suitable framework for such more detailed studies with staggered fermions. 
Eventually one would like to confirm the studies of the equation of state also within 
another discretization scheme, in general with Wilson fermions. Some work in this 
direction is underway [^. However, also here one will have to look for suitable 
improvements of the action in order to reduce the cut-off dependence of thermody- 
namic observables. 
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